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We present a detailed tlieoretical analysis of micro-motion in a time-averaged orbiting potential 
trap. Our treatment is based on the Gross-Pitaevskii equation, with the full time dependent be- 
haviour of the trap systematically approximated to reduce the trapping potential to its dominant 
terms. We show that within some well specified approximations, the dynamic trap has solitary-wave 
solutions, and we identify a moving frame of reference which provides the most natural description 
of the system. In that frame eigenstates of the time-averaged orbiting potential trap can be found, 
all of which must be solitary-wave solutions with identical, circular centre of mass motion in the lab 
frame. The validity regime for our treatment is carefully defined, and is shown to be satisfied by 
existing experimental systems. 



INTRODUCTION 



The Time-averaged Orbiting Potential (TOP) trap 1] 
was an important tool in the realization of Bose-Einstein 
condensates, and it remains a common method for mag- 
netically trapping atoms. Early theoretical descriptions 
of the TOP trap used two approximations: the adiabatic 
approximation, which assumes that the magnetic dipoles 
of the atoms align instantaneously to the magnetic field, 
and the time-average approximation, where the time dy- 
namics of the trapping fields are neglected on the time 
scale of the motion of the trapped atoms. Under these as- 
sumptions, the TOP trap is represented by a static, har- 
monic potential and the condensate eigenstates are rela- 
tively easily calculated (usually by numerical means) and 
are stationary in space. However, condensates formed in 
a TOP trap undergo a spatial micro-motion 0, Q due to 
the underlying dynamic nature of the TOP trap. This 
phenomenon has been studied theoretically under vari- 
ous levels of approximation, by partially lifting the time- 
average approximation Ijj5( or by not applying the adi- 
abatic approximation j^, Q ■ 

In this work we provide a detailed theoretical descrip- 
tion of condensate micro-motion in terms of TOP trap 
eigenstates, including condensate nonlincarity. Our ap- 
proach applies the adiabatic approximation, but partially 
lifts the time-average approximation. Under these condi- 
tions, the TOP trap potential retains some time depen- 
dence and eigenstates of that potential cannot be found 
in the lab frame. However, system eigenstates do ex- 
ist because a frame can be found in which the Hamilto- 
nian for the system becomes time independent. We have 
termed eigenstates of the system found in such a frame 
dynamical eigenstates, since these states are not station- 



ary states in the lab frame. By calculating the dynamical 
eigenstates of the TOP trap, full characterisation of con- 
densate micro-motion is possible. This is essential for an 
understanding of condensate growth and is also required 
for a description of velocity sensitive phenomena occur- 
ing in TOP traps, such as observed in Bragg scattering 
experiments [ 
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In this paper we calculate dynamical eigenstates of the 
TOP trap potential in the quadratic average approxima- 
tion. Within that approximation, the solutions are exact 
in both the linear and the nonlinear case. We begin, in 
section |5J by introducing the TOP trap potential, and 
various approximate forms of that potential. In scction|31 
we derive the transformation to the circularly translat- 
ing frame which we find to be the most natural frame in 
which to investigate the system. In section0|we calculate 
solitary-wave solutions in the quadratic average approx- 
imation and show that dynamical eigenstates calculated 
using the circularly translating frame are a particular 
class of solitary-wave solutions in the lab frame. Iden- 
tifying the dynamical eigenstates of the TOP trap allows 
us to characterize micro-motion and specify the ground 
state of the system. In section [SJ we assess the valid- 
ity of the quadratic average approximation and demon- 
strate that for the typical parameter regime of the TOP 
trap, solitary-wave dynamical eigenstates provide accu- 
rate approximations to the dynamical eigenstates of the 
full TOP trap potential. In section El we discuss lab 
frame solitary- wave solutions which are eigenstates of the 
TOP trap potential in the more commonly used rotating 
frame, and show that these are only a subset of the dy- 
namical eigenstates found using the circularly translating 
frame. We conclude in section [7| 



APPROXIMATE FORMS OF THE TOP 
TRAP POTENTIAL 



number of atoms N, and the s-wave scattering length a, 
i.e. 



The TOP trap consists of a magnetic quadrupole trap 
[i3,ll3 translated by a uniform bias field, whose direction 
rotates at frequency fi Q. For simphcity our discussion 
is presented in a set of dimensionless units defined by the 
position scale xq = y^h/2'muj,j. (a characteristic harmonic 
oscillator length), and the time scale of the inverse of 
the time-averaged trap frequency tUx^ defined below. A 
key feature of the TOP trap is that the zero of the mag- 
netic field follows a circular trajectory of radius rg, and 
trapped atoms are confined well within this trajectory 
(the so-called 'circle of death'), thereby reducing atom 
loss due to spin flips. Typically vq ^ 1000 — 1300 and 

r2~70-i50 0,1113. 



2.1. The adiabatic approximation 

The TOP trap potential in the adiabatic approxima- 
tion is given by 



Vtop ir,t) = rl 



1 



2 [x cos Vtt + y sin Qt) 



To 



y2 + A7J 



(1) 



where r = {x,y,z). That approximation is valid when 
the bias field rotation frequency, f2, is much smaller than 
the Lamor precession frequency [J . 



2.2. The truncated TOP trap potential 



c 



AnhaN 
muJxXQ 



(5) 



2.3. The time-average approximation 

The most common treatment of condensate evolution 
in a TOP trap has also invoked the time-average approxi- 
mation, whereby the potential of equation (0) is averaged 
over a period of the bias field rotation. This leads to the 
time-averaged, truncated form of the TOP trap poten- 
tial. 



VH{r)=rl + ^ix^+y^ + 8z^ 



(6) 



which is a static, harmonic potential, with frequency iOx 
in the x-y plane (in SI units). In the time-average ap- 
proximation, the potential in the time evolution operator 
of equation I^J , is replaced by the time independent trap 
Vh(t), of equation ©. This allows energy eigenstates of 
the system to be readily calculated. 

The time-average approximation is normally assumed 
to be valid when the bias field rotation frequency is much 
larger than the frequency of the time-averaged harmonic 
trap, i.e. in our dimensionless units fi ^ 1. The time- 
averaged treatment neglects system dynamics occurring 
on the fast time scale of the bias field rotation, and it is 
this non-stationary behaviour of a condensate in a TOP 
trap that we describe in this work. 



Expanding the square root of equation (^ in a Taylor 
series, and neglecting terms above second order in the 
small parameter Xa/ro, where Xa is one of x, y, or z, 
leads to the truncated TOP trap potential 



V{r,t) = rl +ro{x cos nt + y sin nt) + -{x'^ +y'^ + 4z^) 



-{xcosV,t + ysinflt)^ 



(2) 



The evolution of the condensate wave function, ■0(r, t), 
is governed by the Gross-Pitaevskii equation. 



i—i>{r,t)^C{r,t)^{r,t), 



(3) 



where for a TOP trap, the time evolution operator, 

£{r, t) = -V2 + V{r, t) + C| V(r, t)\\ (4) 

is time dependent. The truncated TOP trap potential 
F(r,i) is given by equation Q, and C is the dimen- 
sionless nonlinearity coefficient defined in terms of the 



2.4. The quadratic average approximation 

Miiller et al. |2] experimentally observed the dynamic 
effects of the TOP trap on condensate evolution, i.e. 
micro-motion in a TOP trap. Their approach for calcu- 
lating the condensate micro-motion amplitude involved 
balancing the restoring force of the time dependent terms 
of equation |2Jl that are linear in cos Qt or sin 17i, with the 
centrifugal force. In line with this treatment, our work 
invokes what we shall refer to as the quadratic average ap- 
proximation^ where only the terms of equation Q that 
are quadratic in cos Hi or sin U,t are time-averaged. In 
that approximation the TOP trap potential is given by 

Vap(r,t) = VB_{r)+rQ{xcosVLt + ysinnt). (7) 

In the present paper, we calculate dynamical eigen- 
states of the TOP trap potential in the quadratic aver- 
age approximation, where the trapping potential is given 
by equation Q . The accuracy of the quadratic average 
approximation is addressed in section |S1 




FIG. 1: The circularly translating frame, defined by co- 
ordinates R = {X, Y, Z) and equation ^. 



3. THE CIRCULARLY TRANSLATING FRAME 

In the lab frame, the TOP trap potential in the 
quadratic average approximation, given by equation |(7J), 
is time dependent and eigenstates of the Gross-Pitaevskii 
equation can not be found. By transforming to a frame 
that translates in a circular trajectory with radius 7 
(whose value is to be determined) and with angular fre- 
quency fl about the origin of the lab frame, we can re- 
move this time dependence. We refer to that frame as 
the circularly translating frame and we shall see that it 
is the natural frame in which to describe the TOP trap 
system. 

The translation in co-ordinate space is defined by 



R = r — 'J {cos fit, smil.t, 0), 



(8) 



as illustrated in figure^ The momentum in the circularly 
translating frame is derived by differentiating equation 
(|H1), yielding 

P = p+ -7fi(sinf7f, -cosf7t,0), (9) 

where we have used the fact that in our dimensionless 
units p = v/2. 



3.1. Unitary transformation to the circularly 
translating frame 

We now derive the quantum mechanical transforma- 
tion to the circularly translating frame. For clarity, we 
shall denote quantum mechanical operators by O, and 
begin with the linear case. 

The Schrodingcr equation for a single particle state in 
the TOP trap, in the quadratic average approximation, 
is given in the lab frame by 






where 



H. 



p^ + V^p{r,t). 



(10) 



(11) 



The transformation to the circularly translating frame is 
achieved by the unitary transformation 



Uit) = Up{bit))Ur{ait)), 



where 



Uria{t))^e 



- p^P-a{t) 



(12) 



(13) 



translates position by a{t), iJr{a{t))\r) = \r~a{t)), and 

(14) 



Up{b{t)) = e-^^-^(*) 



translates momentum by b(i), ijp{b{t))\p) = |p — b{t)). 
In a comparison with equations © and we find that 



and 



a{t) = 7(cos Q,t, sin Vtt, 0), 



b{t) = --7f7(sinf7i, -cosr2i,0). 



(15) 



(16) 



In the transformation to the circularly translating frame 
lJr{a{t)) and IJp{b{t)) commute, since a{t) ■ b{t) = 0. 
Defining the transformed state vector to be 



W^u{tM), 

equation H10() becomes 



dt 



= H 



ap I 



where 



iJ*p = Up{bit))Uria{t))H,pUl{ait))ul{bit)) 



ip+m) 



da{t) 



db{t) 



(17) 
(18) 

(19) 



dt dt 

The Schrodinger equation in the co-ordinate representa- 
tion can be determined by projecting equation H18(l onto 
state \R). Using the following identities. 



{R\i;)=i;{R,t), 


(20) 


and 




{R\p\i,) = -iyRi>iR,t), 


(21) 


where we have denoted 




/ d d d \ 
^ \dX'' dV dz) ' 


(22) 


yields 




^''^lf'^ - Hl^iR^tmR,t), 


(23) 



where 



da{t) 



R 



zV^ - b{t)) 
db{t) 



dt 



dt 



(24) 



The wave functions in the lab frame and the circularly 
translating frame are related by 



iP\R, t) = e-'^-^'^'^iR + a{t),t). 



(25) 



3.2. Application to the Gross-Pitaevskii equation 

The above derivation, for the Schrodinger equation, 
may also be adapted to the Gross-Pitaevskii equation, 
since the nonhnear term transforms simply under substi- 
tution of equation (|25() . Thus, by substituting a(t) and 
b{t) from equations (|15|l and (|16|l into the Hamiltonian of 
equation H24|l . and including the nonlinear term (which 
is described in terms of the new density |V'*(-R, i)P), the 
Gross-Pitaevskii equation in the circularly translating 
frame is 



where 



and 



i^^i:'{R,t) ^ Cl^{R,t)i,\R,t), 



CUR,t)^HlJR,t) + C\i;\R,tf 



(26) 



(27) 



Hl{R,t) = -V^ + FH(i?) 



- (7 + 2ro - 7f^^) (X cos nt + Y sin nt) 



+7ro--7'(f^'~l). 



(28) 



The single particle Hamiltonian of equation (|28|l is iden- 
tical to the Hamiltonian derived using a classical frame 
transformation to a non-inertial frame of reference [l5j, 
applied to the circularly translating frame. Choosing 
7 = 7t, where 



7t = 



2ro 



the evolution operator of equation H27|) simplifies to 



(29) 



4p(«,0 



R 



ViiiR) + C\ij\R,t)\^+e, (30) 



where 



e=-j^{n^-l). 



(31) 



The energy offset e can be interpreted by expressing 
equation l|31(l in the form 



yap([7t,0,0],t-0)-lij(O)-Sf, 



(32) 



The first two terms represent the additional potential 
energy due to the displacement of a point body from 
the trap centre to radius 7t. The remaining term, E^i ~ 
7^0^/4, represents the energy of a point body rotating 
about the origin of the lab frame, at a radius 7t with fre- 
quency 57, which is simply the expected energy shift as- 
sociated with the transformation to the circularly trans- 
lating frame |15l |. 

The time evolution operator in the circularly translat- 
ing frame, as given by equation (|30|l . substituted into 
equation H2()|) yields the Gross-Pitaevskii equation for a 



time independent harmonic trap, with an energy offset 
of e. Thus, eigenstates of the TOP trap in the circularly 
translating frame exist in the quadratic average approx- 
imation. For clarity we write these as 



xP\R,t) ^ iPiiiR)e 



-i(fiii+e)t 



(33) 



where 'ipii{R) are the well-known solutions to the time 
independent Gross-Pitaevskii equation for a time inde- 
pendent harmonic trap, i.e. 



^hV'h(-R) 



R 



Vn{R) + C\iJn{R)\^ MR) 



(34) 



3.3. Generalization to quantum field theory 

The transformation given by equation (|25|l can be ap- 
plied to the operator Heisenberg equations of motion for 
the full quantum field operator "^{r^t). In much the same 
way as our discussion above, this yields the equation of 
motion in the circularly translating frame 



. d-9tiR,t) 
'' dt 



-V|^ + Vn{R) + Ci'l{R,t)i'tiR, t) 
-e]^t(H,t). (35) 



Since this represents the full quantum field theory, the 
motion of uncondensed particles is also correctly treated 
in the circularly translating frame. 



4. SOLITARY- WAVE SOLUTIONS 

Solitary- wave solutions, where the wave function 
evolves without changing shape, can be found for the 
TOP trap in the quadratic average approximation. Mor- 
gan et al. [1^ have shown that the Gross-Pitaevskii equa- 
tion, with particular forms of potential, has solitary- wave 
solutions which propagate in one dimension of a multi- 
dimensional space. That work was extended by Mar- 
getis ^13j where solitary-wave solutions may have center 
of mass motion in any of the space dimensions. Also, 
Japha and Band 14] have shown that in a moving har- 
monic trap the motion of the condensate centre of mass 
can be entirely decoupled from the evolution of the con- 
densate shape. We have extended the derivation by Mor- 
gan et al. [ij to include the case where solitary-wave 
solutions can propagate in three dimensions, as was indi- 
cated to be possible by Margetis [Tjl . In the following we 
present a brief summary of the results of our derivation. 

We begin by postulating that solitary-wave solutions 
to the TOP trap will have the form 



^sw(r,t) = ^H(r--f(i))e-*^H*+^^(^'*), 



(36) 



where the envelope wave function ^h(^) is an eigenstate 
of the time independent Gross-Pitaevskii equation for the 



TOP trap potential in the time-average approximation 
with chemical potential /in, i-e. defined by equation (|34|l . 
The position offset in the envelope wave function is 



f{t) 



ilj^{r)ril;}i{r)dr, 



(37) 



which can be interpreted as the time dependent position 
of the centre of mass of the solitary- wave since the second 
integral is zero due to the particular form of Vh(t"). The 
phase S{r^ t) is determined by substituting the solitary- 
wave solution, equation (|36(l . into the time dependent 
Gross-Pitaevskii equation, equation (Q, where C{r,t) is 
replaced by Cap{r, t), in which the quadratic average ap- 
proximation is used, i.e. 



'Cap(r,i) 



Kp(r,i) + C|V(r,i)|2 



(38) 



where T4p(r,t) is given by equation Q. Taking a sim- 
ilar approach to that of Morgan et al. [l2| the Gross- 
Pitaevskii equation can be separated into real and imag- 
inary parts yielding two equations. The equation derived 
from the imaginary part can be simplified by writing 



, ^ 1 df(t) , , 



(39) 



We choose the trivial solution, K{r,t) — K{t), which 
is the only possible solution in the one dimensional case 
[l^ and has also been suggested as the unique solution 
in the general case ^3- Substituting the trivial solution 
into the equation derived from the real part of the Gross- 
Pitaevskii equation we find that 



m^ll 



'(f)- 



m" 



dt. 



(40) 



By equating mixed differentials of S{r,t), the centre of 
mass motion of the solitary-wave solutions can be found 
to obey 



1 d'^f{t) 

2 dt^ 



-VF(r,i), 



which is a form of Ehrenfest's theorem, where 



F(r,i) = Kp(r,t)-FH(r-f(t)). 



(41) 



(42) 



For solitary- wave solutions to exist, in the form that we 
have discussed, both sides of equation (|41|) must be in- 
dependent of r and therefore the function F{r,t) must 
be at most linear in r. The TOP trap potential in the 
quadratic average approximation obeys this criteria and 
thus solitary-wave solutions exist. It is possible to solve 
for f{t) which has six constants of integration, given 
by the initial values of the centre of mass position and 
momentum of the particular solitary-wave solution (see 
equation IjAip l. 



4.1. Solitary-wave solutions which are eigenstates 
in the circularly translating frame 

The dynamical eigenstates calculated using the circu- 
larly translating frame, given by equation H33|l . are a par- 
ticular class of solitary-wave solutions in the lab frame. 
This can be confirmed by transforming the solitary-wave 
solutions, as given by equation (|36|l . into the circularly 
translating frame, and requiring that these solutions sat- 
isfy the time independent Gross-Pitaevskii equation in 
that frame, i.e. 



where 



and 



^^lw^swiR^t) ^ ^LiR,mlwiR^t) 



^'swiR.t) = VH(r-r(i))e-''^H*+'^(^'*) 



e2 



Msw 



ii7tr2r-(sinnt,-cosr2t,0) 



Mh + £• 



(43) 



(44) 



(45) 



In order to satisfy equation H43|l . the solitary- wave solu- 
tions of equation (|44|l have a restriction on f(i), as de- 
rived in appendix ^ Solitary- wave solutions of the lab 
frame are eigenstates of the TOP trap in the circularly 
translating frame if and only if the initial conditions of 
the centre of mass motion of the solitary-wave solutions 
have particular values such that 



r{t) = 7t(cosr2t, sinfJijO), 



and therefore 



Pit) 



1 



7trj(sinrJi, — cosftt, 0). 



(46) 



(47) 



These equations represent circular motion at the TOP 
trap frequency fl, with radius 7t. Previously 01 13 ^^^ 
micro-motion position amplitude has been determined to 
be 2ro/fi~, which is in agreement with our result (see 
equation H29|) ) in the limit fl ^ 1. 

With the centre of mass motion for solitary-wave dy- 
namical eigenstates of the TOP trap given by equations 
(gBJ and 1(171) . we find that R = r - r{t) and the phase 
S{r, t) becomes 



S{r,t) 



1 



7tf}r • (sinm, -cosrit,0) -rf. (48) 



Making these substitutions, equation (|44|l simplifies to 
equation H33|l so that all dynamical eigenstates of the 
TOP trap, calculated using the circularly translating 
frame, are a particular class of solitary-wave solutions in 
the lab frame with centre of mass motion given by equa- 
tions H46|l and H47I) . This shows that the origin of the 
circularly translating frame (refer to equations © and 
@ with 7 = 7t) moves with the centre of mass motion of 
the solitary- wave dynamical eigenstates of the TOP trap, 
therefore justifying our choice of the circularly translat- 
ing frame for describing the TOP trap. 




10 



y 



FIG. 2: A particular eigenstate of the two dimensional equiv- 
alent of equation 13411 . calculated numerically. The two di- 
mensional nonlinear strength is C2D ~ 600 and the chemical 
potential is hh = 10.56. 



4.2. Dynamical eigenstates 

All the dynamical eigenstates of the TOP trap follow 
the same circular trajectory in the lab frame, as given by 
equations l|^ and l|Tf|l . This motion is independent of 
both the chemical potential of the state and the nonlinear 
strength of the system. Furthermore, the solitary-wave 
dynamical eigenstates retain their orientation with the 
lab frame throughout their trajectory. 

As an example, a two dimensional excited state of the 
envelope wave function ipn^x, y), with a nodal line along 
the y axis, is presented in figure |2| In the lab frame, the 
solitary-wave dynamical eigenstate of the TOP trap, cor- 
responding to the envelope wave function in figure 1^)1. 
consists of the envelope wave function moving in a circu- 
lar trajectory, while maintaining its orientation with the 
lab frame, and the orientation of the nodal line along the 
y axis. 

In the linear case (C — 0) the solitary-wave dynamical 
eigenstates form a complete basis for the TOP trap in the 
circularly translating frame. This is because these states, 
in the circularly translating frame, are eigenstates of the 
harmonic oscillator equation with an additional energy 
offset (see equations (PUJI and (|23))- 



malize during evaporation into the usual Bose-Einstein 
distribution, and hence condensation from the vapour 
will be into the TOP trap ground state, as determined us- 
ing the circularly translating frame. Therefore, the state 
into which bosons condense, in the quadratic average 
approximation is the solitary-wave dynamical eigenstate 
given by equation l|44(l with the envelope wave function 
being the ground state of equation H34|l and f{t) given 
by equation (^BJ. 



5. VALIDITY OF APPROXIMATIONS 

5.1. Corrections to tiie quadratic average 
approximation 

Throughout this work we have used the quadratic av- 
erage approximation without assessing its validity. Here 
we give a systematic assessment of the validity regime of 
the quadratic average approximation for the linear case 
(where the mean- field interaction is neglected). This al- 
lows simple analytic results to be obtained. 

The single particle Hamiltonian, with the truncated 
time dependent TOP trap potential of equation ((2Jl, takes 
the form (in the circularly translating frame) 



H' = Hl^ + WiR), 



(49) 



where iJ'„ is the single particle Hamiltonian in the 



ap 



quadratic average approximation, i.e. 



HL = P 



ap 



X 



P^ + Pl + Vu{R)+e. 



(50) 



In equations (|49() and (|50|l we have used an operator for- 
malism where the position and momentum component 
operators are denoted by R and Pj, respectively. The 
perturbative potential, W{R), accounts for the remain- 
ing terms of the TOP trap potential of equation ^ that 
are not retained in the quadratic average approximation. 
In the circularly translated frame, these terms are given 
by 



W{R) 



-{X^ - Y^){cos^ nt - sin^ nt) 
-7t {X cos nt + Y sin fit) 



4.3. Condensation and the ground state of the 
TOP trap 

The ground state of the TOP trap system is the 
solitary-wave dynamical eigenstate with the lowest chem- 
ical potential in the circularly translating frame. This 
occurs when /^h, the chemical potential of the envelope 
wave function, takes its lowest possible value (see equa- 
tion H45|) '). Since, as noted in section rOl the uncon- 
densed atoms experience the same potential as the con- 
densate, in the circularly translating frame, these ther- 



-XY sin ilt cos ilt 7^ . 

4 "■ 



(51) 



The harmonic oscillator creation and annihilation op- 
erators in the circularly translating frame, are defined in 
our dimensionless units, as 



1 



ax = -X + iPx 


(52) 


ay = 2^ + iPY 


(53) 


az = V2Z + iPz, 


(54) 



where [aj, a]^] = Sjk, and J and K are one of X, Y, or 
Z. Making these substitutions we find that 



Hi 



and 



ax<^x + al^ay + a\az + e + Tq + 1 + \2, (55) 



,t2 



W{R) = -^(a^^ + a^ + 2a5(^ajf)(cos^f7i-sin^rJi) 
+ -(a|? + fly + 2a|,aY)(cos^ fit - sin^ ilt) 
— -7t(ax + ajc) cosVtt — -7t(ay + ay) sinfJi 



(a^ + ax ) (oy + fly ) sin fit cos J7i 



1 2 



(56) 



Utihsing the number operator kets, which satisfy 
ajaj\nj) — nj\nj)^ the eigenket of the single particle 
Hamiltonian is |nx7?^F7?^z)', i-e. 

fl'*p|nx,ny,nz)' ^ E^\nx,nY,nzY . (57) 

The energy spectrum is given by 

E^ ^nx+nY + nz+e + rl + l + V2, (58) 

in agreement with equation H45|l . The energy spectrum 
of the harmonic oscillator terminates at the ground state, 



|0, 0, 0)', which has energy £:*=£■ 



V2. 



Using time dependent perturbation theory, the evolu- 
tion of the ground state to first order in the perturbation 
W{R) is given by 



W = A{t). 



,-'(4-i7t)t 



[|0,0,0)* 



piQ,t p — iQt 



|1,0,0)' 



7t / 2e-'^ 

'4 \n^ - 1 ^ n + 1 n-1 

1 / 2e-2«* 



2if2t -2iat 



-|0,2,0)*)], (59) 

where A{t) is a constant of normalisation. From this 
expression we can deduce that the quadratic average ap- 
proximation is valid in the linear case, within the param- 
eter regime where 7t ^ 51 and 1 <C il. 

5.1.1. Nonlinear case 

It is clear that a perturbative two-timescale asymptotic 
expansion in powers of l/Jl could be made for the non- 
linear case. Thus, we expect that for the nonlinear case, 



the quadratic average approximation is also valid within 
the regime derived above for the linear case. 

We have carried out two dimensional numerial calcu- 
lations for the nonlinear case which verify this. For ex- 
ample, for a typical TOP trap system where tq = 1241 
and Vt ~ 153, we have propagated the Gross-Pitaevskii 
equation using both the truncated TOP trap potential 
of equation |2Jl and the potential in the quadratic av- 
erage approximation, given by equation |7J. An appro- 
priate value for the two dimensional nonlinear stregth 
is C2D = 600 which correspondes to iV ~ 2x10^*^ in the 
Otago TOP trap [31 ■ The initial state was chosen to be 
the ground state of the TOP trap in the quadratic aver- 
age approximation, i.e. the ground state of equation l|M|) 
calculated numerically using optimization methods and 
shifted in position and momentum according to equations 
(jIH|l and (|T7|) (at t — Q). That state was propagated 
by the Gross-Pitaevskii equation for one period of the 
bias field rotation for two cases: (i) with the truncated 
TOP trap potential, giving V''™'"'(a:,j/, i = 27r/17), and 
(ii) with the potential in the quadratic average approxi- 
mation, yielding tp'''^'^'^ {x , y,t = 2-K/Vt). The method used 
was an accurate fourth order algorithm, with a grid of 
512x512 points over a 60x60 range in position, and 20000 
time steps. The deviation between the two solutions was 
found to be / |V'*™"'^(a;, y, t = 2-K/n) - iP'i'''"^ {x , y , t = 
2t: / n)\^ dxdy = 4:.44:xl0-^. 



5.2. Validity of solitary-wave dynamical eigenstates 

The validity of our solitary- wave dynamical eigenstates 
as dynamical eigenstates of the full TOP trap depends on 
three validity conditions: 

1. the adiabatic approximation, 

2. the truncation of the TOP trap potential from 
equation ^ to equation Q, and 

3. the quadratic average approximation. 

The validity of the quadratic average approximation was 
addressed above and applies in the regime where 7t <C il 
and 1 ^ ri. We note that the condition jt <^ il can 
be re- written as 2ro ^ fi'^. The truncation of the TOP 
trap potential to yield V{r,t) of equation |(2J) is valid 
provided Xa <C tq where Xa is one of a;, y, or z. A useful 
estimate of Xa is given by the sum of the Thomas-Fermi 
radius of the solitary- wave dynamical eigenstate envelope 
wave function with the lowest chemical potential, and the 
radius of the dynamical eigenstates trajectory, "ft. This 
yields 



/3ocy 

\V2^) 



■7t, 



(60) 



where C is given by equation ^ . The adiabatic approx- 
imation is valid when the bias field rotation frequency fl 
is much smaller than the Larmor precession frequency. 
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given in our dimensionless units by the potential. As an 
estimate of the Larmor precession frequency we use the 
magnitude of Vtop('";0 which, assuming that Xq <C tq, 
is of the order of r^. Thus, the adiabatic approximation 
is vahd provided that il ^ r^. 

Finally, collating the validity regimes we find that our 
solitary-wave dynamical eigenstates calculated using the 
circularly translating frame, are an accurate description 
of the dynamical eigenstates of the full TOP trap system 
within the parameter regime where 



1 < n 

30C < \/27rr^. 



(61) 



Typical experimental parameters are well within these 
criteria. As an example, the Otago TOP trap system of 



Rubidium-87, where a — 55x10 



and 



= 18 Hz 



[9j, leads to the third validity condition from equation 
i|()T|l becoming N <^ Tq so all three conditions are easily 
satisfied. 



6. THE ROTATING FRAME 

In previous theoretical work the rotating frame has 
been used to calculate eigenstates of the TOP trap sys- 
tem under various levels of approximation [J, |a| • How- 
ever, it can be shown that using the rotating frame to 
describe the TOP trap allows only a limited set of dy- 
namical eigenstates to be found. For completeness, we 
present this calculation in appendix IbI where, using the 
same methods as in sections|31and0I we show that a par- 
ticular class of solitary-wave solutions in the lab frame 
are eigenstates of the time independent Gross-Pitaevskii 
equation in the rotating frame (equation (|B6p l. As be- 
fore, these solitary-wave dynamical eigenstates follow a 
circular trajectory in the lab frame, described by equa- 
tions (|46() and H47() . but unlike solitary- wave dynami- 
cal eigenstates calculated using the circularly translating 
frame, the solitary-wave solutions which are eigenstates 
of equation (jB6|l must also obey an additional symme- 
try, which is that the envelope wave function must be an 
eigenstate of the z component of angular momentum (see 
equation (|B8I) ). This condition requires solitary- wave dy- 
namical eigenstates calculated using the rotating frame 
to be cylindrically symmetric about their centre of mass. 
Figure |21 shows an example of a dynamical eigenstate 
envelope wave function with a corresponding solitary- 
wave dynamical eigenstate which does not satisfy the 
time independent Gross-Pitaevskii equation in the rotat- 
ing frame. Physically we can see why: the nodal line of 
that solitary-wave dynamical eigenstate, which remains 
oriented along the y axis in the lab frame, will appear 
to rotate in the rotating frame so that the dynamical 
eigenstate is not stationary in that frame. 

Solitary-wave dynamical eigenstates of the TOP trap 
retain their orientation with respect to the lab frame as 



they move. Consequently, the rotating frame is not an 
appropriate choice for describing dynamical eigenstates 
of the TOP trap system, because the eigenvalue equa- 
tion in the rotating frame incorporates the angular mo- 
mentum operator, and places additional symmetry con- 
straints on dynamical eigenstates of the system that are 
not in general necessary. 



7. DISCUSSION 

We have carried out a detailed characterization of con- 
densate micro-motion in a TOP trap, under some well de- 
fined approximations. These approximations, which are 
well justified for typical TOP traps, are (i) the adiabatic 
approximation (which neglects spin precession effects), 
(ii) the assumption that the condensate is located well 
within the circle of death, and (iii) the quadratic average 
approximation (which time-averages quadratically oscil- 
lating terms in the potential). Our treatment allows for 
condensate nonlinearity and we have shown that within 
these approximations, solitary-wave solutions of the non- 
linear Gross-Pitaevskii equation exist. We have identified 
the circularly translating frame as the most appropriate 
frame for describing the system, and have shown that 
eigenstates can be found in that frame, and that they 
must all be solitary-wave solutions of a certain type. In 
particular, all of the solitary-wave dynamical eigenstates 
have identical centre of mass motion, which in the lab 
frame is a circular trajectory with radius "ft and momen- 
tum magnitude oTt^- 

Previous theoretical discussions of dynamical eigen- 
states of the TOP trap have been given within similar 
approximations, but with the additional restriction that 
the nonlinearity due to the atomic interactions is either 
approximated or neglected. Kuklov et al. |j| have ob- 
tained exact eigenstates for the linear Schrodinger equa- 
tion within the adiabatic approximation using the trun- 
cated TOP trap potential of equation (|2Jl and those au- 
thors have also presented an approximate many-body 
treatment. Their exact single particle solutions are ob- 
tained using numerous transformations and the form of 
eigenstate micro- motion is not readily evident. Their 
method also employs the rotating frame which, as we 
have shown within the quadratic average approximation, 
limits the possible dyiiamical eigenstates that can be 
found. Minogin et al. [a] have used an approximate inter- 
action picture method which provides information about 
the atomic momentum modulation in a TOP trap, but 
does not describe the micro-motion in the position co- 
ordinates. 

Our choice of the circularly translating frame allows 
solitary-wave dynamical eigenstates, which retain their 
orientation relative to the lab frame, to be readily iden- 
tified for the TOP trap system. These dynamical eigen- 
states have no restriction on the z component of angu- 
lar momentum of the envelope wave function. By con- 
trast, we have shown that the dynamical eigenstates cal- 



culated using the rotating frame constitute only a subset 
of the dynamical eigenstates calculated using the circu- 
larly translating frame, and are required to be cylindri- 
cally symmetric about their centre of mass. 

Finally, we have shown that the validity regime for 
the quadratic average approximation is defined by the 
conditions 1 <C ^ and 2ro <C ^^- These criteria are well 
satisfied by existing TOP trap systems. 
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APPENDIX A: RESTRICTIONS ON 

SOLITARY- WAVE SOLUTIONS WHICH ARE 

EIGENSTATES IN THE CIRCULARLY 

TRANSLATING FRAME 

We seek solutions to equation 143|l which have the form 
of equation (|44|l . We will derive the particular form for 
the centre of mass motion, f{t), for which such solutions 
exist. The centre of mass motion can be solved in general 
using equation (|41|) . which gives 

x{t) — (xi — 7t) cosi -|- vi sini -I- jt cos fit 
y{t) = X2 cost + {v2 — 7tfi) sint -I- jt sinflt (Al) 
z{t) = X3 cos 2V2t + -^ sin 2\/2i, 



2%/2 



{xi,X2,X3) and 



where we have defined r\u^Q\ 

df{t)/dt\(t=Q) = {vi,V2,Vi). 

Eigenstates of the circularly translating frame must 
satisfy both equation (|43() and 



d^l^{R,t) _ 



dt 



/^sw^sw(-R>i)- 



(A2) 



Substituting the solitary wave solution of equation 144() 
into equations (|43|l and (|A2p . and making the change of 
variables s = r — f{t), yields respectively 



1 



1 
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Msw = m + T^it^' + -y{t) + -y{t) + 



df{t) 



dt 



--7tf(t) • (cosm, sinf^i, 0) -f -s ■ r{t) 



2 ' \ dt 



(sin fit, — cos ilt, 0) 



7 1 
+ -S ■ (0, 0, z{t)) 7ts • (cos nt, sin nt, 0) 



df{t) 



dt 



7til(sin fit, — cos fit, 0) 



9(«)(A3) 



and 



A'sw 



A^H + ^7t ^' - ir'it) 



2 ^^^ dt^ 



1 fdf{t) 
4 V dt 



- --/t^'^r{t) ■ (cos m, sin f7i, 0) 



4-K^) 



(sin fit, — cosftt, 0) 



1 
'2^ 



dt^ 
df{t) 



-^tVt^s ■ (cos J7i, sin fit, 0) 



dt 



^tVl{s\\\Vlt, — cosflt, 0) 



where 



g{s) 



VsV'H(g) 

V'h(s) 



■9{slkA) 



(A5) 



Equating equations l|A3p and l|A4p . and substituting 
the general form of f{t), given by equation (jAljI . we find 
that equality requires xa = and W3 = and therefore 
z{t) = 0. Inserting this result, and the general expressions 
for x{t) and y[t) from equation HA1|) . equations IjASp and 

^41) both simplify to 



Msw - Mh 
1 



i7?o' 



+ j{xl +xl+vl+vl- 27t(a;i + V2fl)) 
+ig{s) ■ ((.Ti — 7t) sini — vi cost, 
X2 sin t — {v2 — 7t^) cos t, 0) 



-2*- ((^1 



7t) cosi + fi smi, 
X2Cosi + {v2 — 7t^) sini,0). 



(A6) 



For the solitary-wave solutions of equation H44|l to be 
eigenstates of the TOP trap in the circularly translating 
frame, the chemical potential /i^w must be independent 
of s and t. Solving equation 
t = tt/2, we find that 



Inipnis) 



1 

--;;iSiS2- 
2 t;ia;2 



[6|l for g{s), at t = and 
Xl + {V2 - 7tf^)^ 



-lSiS2 

2 wia;2 

+^2(53), 



Xl - 7t)(w2 

(xi -7t)^ +vf 



7tfi) 



(Xl 



7t)(w2 - 7tf^) 



+ Ci(s3) 



(A7) 



where s = (si, 32,53), and 6*1(53) and 6*2(53) are con- 
stants of integration. The only possible solution is there- 
fore (xi, X2, wi, ^2) = (7t, 0, 0, 7tr2), which illiminates 
g{s) and s from equation ljA6p yielding 



Msw — A^H 



i7t("^ 



1), 



(A8) 



which is in agreement with equation (|45|l . Concluding 
then, solitary-wave solutions described by equation 1)44(1 
which are eigenstates of the TOP trap in the circularly 
translating frame exist if and only if f\(^t=o) — (7t,0, 0) 
and (if(t)/(ii|(f=o) = (0,7t^, 0), and therefore the centre 
of mass motion of these states in the lab frame, given in 
general by equation I|A1(I . simplifies to 



r{t) — 7t(cosr2i,sinOi,0). 



(A9) 
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APPENDIX B: SOLITARY- WAVE DYNAMICAL 

EIGENSTATES DERIVED USING THE 

ROTATING FRAME. 

In this appendix we calculate dynamical eigenstates 
of the TOP trap potential in the quadratic average ap- 
proximation using the rotating frame. That frame, with 
co-ordinates r' — {x',y',z'), rotates at the frequency of 
the bias field, and is defined by the co-ordinate transfor- 
mation 



X — X cos fit + y sin Qt 
y' — —xsmi^t + ycosftt 



(Bl) 



In the rotating frame the Gross-Pitaevskii equation be- 
comes 



^^V/(r',0 



(B2) 



where the evolution operator in the rotating frame is time 
independent and is given by 



CiJr',t) = C,p{r,t)-nL{r') 



VlJr')~nUr') 



+C\^'{r',t)\' 



(B3) 



and the wave function in the rotating frame is ip'^ir' , t) = 
^(r, t). The angular momentum in the rotating frame has 
a component in the z direction given by 



Ur') 



y 



_d_ 
dx' 



'A 

dy' 



(B4) 



and we note that L^{r') = L,^{r) IT3. Finally, the TOP 
trap potential of equation Q becomes, in the rotating 
frame, 



VUr') 



l{x' + 2ror + \{y'' + 8z'% 



(B5) 



which is a stationary harmonic potential shifted from the 
origin. 

Eigenstates of the TOP trap in the rotating frame obey 
the time independent Gross-Pitaevskii equation in that 
frame. 



fi^Y(.r\t)^C:Jr\t)i;'^ir',t). 



(B6) 



Here we show that a particular class of solitary-wave so- 
lutions obey equation HB6|) . We denote solitary- wave so- 
lutions which are TOP trap eigenstates in the rotating 
frame by '0sw('"'' '')' with chemical potential fJ-sw Trans- 
forming the solitary-wave solution of equation (|36f) into 
the rotating frame yields 

^i„'(cosOt^+sinnti^) 
i,y(cosnt^-sinmi^) 



'(^) 



(B7) 



The detail of substituting equation ljB7|l into equation 
(|ij6|) and enforcing /i"^ to be independent of spatial and 
temporal co-ordinates is given in section ^ of this ap- 
pendix. The results are discussed here. 

Solitary-wave dynamical eigenstates of the TOP trap 
potential, in the quadratic average approximation, as cal- 
culated in the rotating frame, have two restrictions. The 
first is that the centre of mass motion of the dynamical 
eigenstates must be given by equations (jIH|l and l|T7|l . 
This is not suprising since we also found this restriction 
on solitary-wave dynamical eigenstates calculated using 
the circularly translating frame. The second restriction 
is that 



Lz{r)'ipii{r) =l^'ijjii{r), 



(B8) 



enforcing the envelope wave function to be an eigenstate 
of the z component of angular momentum. This means 
that dynamical eigenstates, as calculated using the rotat- 
ing frame, must have a cylindrically symmetric density 
about their centre of mass. We found in our discussion of 
dynamical eigenstates in section !?!^ that non-symmetric 
dynamical eigenstates do exist for the TOP trap and that 
they are solitary- wave solutions in the lab frame retaining 
their orientation to that frame. Solitary-wave dynamical 
eigenstates which also satisfy equation ljB8|) are only a 
subset of the dynamical eigenstates found using the cir- 
cularly translating frame. 

Substituting equations (j^ and HB8|) into the solitary- 
wave solution in the rotating frame, equation ljB7p . the 
solitary-wave dynamical eigenstates calculated using the 
rotating frame, i.e. satisfying equation ljB6|) . have the 
form 



i^l^{r\t) ^ ^{r')e 



= fi^('r''Vi*')'tf^y'-Vsw* 



(B9) 



where we have explicitly extracted the time dependence 
in the envelope wave function by writing 

^Puir - f(i)) = 0(r')e'™'=*. (BIO) 

The chemical potential spectrum in the rotating frame is 

/^sw = /^H +£-f^^z, (Bll) 

where Iz is defined by equation (|B8|I . The wave function 
phase 7tOy'/2, in equation ljB9|) . is derived from the co- 
ordinate dependent phase of equation ljB7|) and accounts 
for the centre of mass momentum of the eigenstates, given 
in the lab frame by equation (|47|l . The chemical potential 
spectrum of the rotating frame, given by equation I)B11|I , 
can be decomposed into three parts: the energy of the 
state that forms the envelope, the additional energy off- 
set £, and an angular momentum term arising from the 
rotating frame, as expected [l5| . 



Restrictions on solitary-wave solutions which are 
eigenstates in the rotating frame 



We seek solutions to equation (|B6p which have the 
form of equation (jB7|) . The conditions required for such 



11 



solitary-wave solutions to be eigenstates of the TOP trap 
in the rotating frame can be found foUowing a similar 
structure to that used in appendix El for the circularly 
translating frame. 

Solitary-wave solutions which are to be TOP trap 
eigenstates in the rotating frame must satisfy both equa- 
tion IIB6II and 






(B12) 



Substituting the solitary wave solution of equation ljB7p 
into equations (|B6I) and ljB12p . and making the change 
of variables s = r — f(t), yields respectively 



Msw = A*H - -^r{t) 



dH{t) 
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+f{t) ■ roicosnt,smnt,0) + -s ■ (0,0, z(i)) 
-I — s ■ [f{t) + 2rQ{cosnt,smnt,0)] 



2 \ dt J 

2 \ dt dt 



and 

Msw 



A^H 



dfjt) 
dt 



-Mt) 



n{y{t),-x{t),0) 



4 
g{s), (B13) 



(ff(t)_ 
~d£^ 



:-r\t) 



1 d'^rjt) 
2^' dt^ 



4a(^),«„,-.„,„,.l(^)- 



2 \ dt dt 



niJs) 



df{t) 



dt 



17(y(t),-5;(t),0) 



9is), (B14) 



where g{s) is given by equation IjASp and 



Us) 



■0h(s) 



(B15) 



Equating equations l|B13p and l|B14|l . and substituting 
the general form of f{t), given by equation HA1|I . we find 



that equality requires 2:3 = and f 3 = and therefore 
z{t) — 0. Inserting this result, and the general expressions 
for x{t) and y{t) from equation (|A1|I . equations l)B13|l and 
(|B14p both simphfy to 

Msw = m - ^lz{s) + -{xl + xj + vl + vi) 

+ -{fl{viX2 - W2a;i) + ItXiifl'^ - 1)) 

—ig{s) ■ (ci cost + C2 sint, C3 cosi — C4 sini, 0) 

1 , 
^t:^ ■ (C2 cost — ci smt, 

—C4 cost — C3 sin t, 0), (B16) 

where ci — vi + nx2, C2 — V2fl — xi — 7t(^^ — 1), 
C3 = f 2 — ^xi, and C4 — vifl -\- X2- For the solitary- wave 
solutions in the lab frame to be eigenstates of the TOP 
trap in the rotating frame, the chemical potential /ig^ 
must be independent of s and t. Solving for g{s), at t = 
and t = Ti 12, we find that the only possible solution oc- 
curs when c\ = C2 and C3 = — C4. Substituting back into 
equation (|B16I) we find that (01,02,03,04) = (0,0,0,0), 
or rather {x\,X2,v\,V2) = (7t, 0, 0, 71^2). Furthermore, 
we find that ^2(5) must be independent of s, and the 
chemical potential in the rotating frame is 



1 



Msw = A^H + -.It (^ - 1) - ^'^z 



(B17) 



which is in agreement with equation (|B11|I . Concluding 
then, solitary- wave solutions which are eigenstates of the 
TOP trap potential in the rotating frame exist if and only 
if r|(t^o) = (7t,0,0), dr(t)/dt|(t=o) = (0,7tf^,0), and 



Lz(s)V'h(s) = ^2'0h(s). 



(B18) 



The centre of mass motion in the lab frame of solitary- 
wave dynamical eigenstates calculated using the rotat- 
ing frame is identical to that of solitary-wave dynamcial 
eigenstates found using the circularly translating frame. 
Equation (|B7p can now be written as 



^l^(r\ t) = Mr- f(t))e5'^*"^'-*(^H+e)*^ (g^g) 



and using equations (jB12|) and (|B18|) it can be shown 
that 



V'H(r--r(t)) = (j){r')e 



l\iQ,l^t 



(B20) 
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